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Abstract 

In quantum field theory (QFT), the vacuum expectation value of a normal product of creation and 
annihilation operators is always zero. This simple property paves the way to the classical treatment 
of perturbative QFT. This is no longer the case in the presence of initial correlations, that is if the 
vacuum is replaced by a general state. As a consequence, the combinatorics of correlated systems 
such as the ones occurring in many-body physics is more complex than that of quantum field theory 
and the general theory has made very slow progress. Similar observations hold in statistical physics 
or quantum probability for the perturbation series arising from the study of non Gaussian measures. 
In this work, an analysis of the Hopf algebraic aspects of quantum field theory is used to derive the 
structure of Green functions in terms of connected and one-particle irreducible Greeen functions for 
perturbative QFT in the presence of initial correlations. 

1 Introduction 

In quantum field theory (QFT), the initial state is most often the vacuum. Many quantum field concepts, 
such as Feynman diagrams, the Dyson equation and the Bethe-Salpeter equation rest on the special 
properties of the vacuum. These desirable concepts extend to special states called quasi-free states 

For general states, it is not possible to write the Green function in terms of standard Feynman 
diagrams and the structure of the Green functions is more complex. For example, the Dyson and Bethe- 
Salpeter equations do not hold. The Dyson equation describes the structure of the two-point interacting 
Green function. Its extension to non-quasi-free states was discovered by Hall [J] in 1975. However, the 
equivalent structural equations for 2n-point interacting Green functions (with n > 1) is not know for 
general initial states. The determination of this structure is the main purpose of the present article, with 
a particular emphasis on the notion of one-particular-irreducibility in this framework. 

We stress that the problem of the calculation of Green functions for initial states that are not quasi- 
free has important applications. For instance, many highly-correlated materials contain transition metals 
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where states of the 3d shell are degenerate. The consequence of this degeneracy is that a small external 
perturbation can create a very strong change in the state of the system. For instance, a small external 
magnetic field induces a strong variation in resistance (giant magnetoresistance), that is used to build 
high-density storage disks. The knowledge of the Green functions would enable us to calculate accurately 
the properties of such materials. 

Although we will use mainly the langage of QFT, let us mention that the situation where the vacuum 
is the initial state corresponds in statistical physics to the case of Euclidean measures as showing up e.g. 
in Euclidean Quantum Field Theory. This is because of the joint use of the Wick theorem (in its simplest, 
Gaussian, form). For general measures, the problem of determining the fine combinatorial structure of 
perturbation series was addressed in [S] where, in particular, the combinatorics of truncated (or connected) 
moment functions was studied intensively. We refer to this article, also for further motivations and 
examples of computations of thermodynamic limits involving the use of generalized Feynman diagrams 
as we also use them. 

Concretely, in the present article, we investigate the structure of Green functions with Hopf algebraic 
methods. Hopf algebras have been implicitly used for a long time in statistical physics and quantum field 
theory: Ruelle 6 , Borchers T and Stora [5] used a product that is called the convolution product in 
the Hopf language. Wightman and Challifour 9J defined a triple dot product that was rediscovered only 
much later in quantum chemistry [lOj and was given an algebraic meaning in [llj . where it was called 
adapted normal product. The work of Wightman and Challifour was summarized by Stora in ref. [8]. 

The explicit introduction of Hopf algebras by Kreimer and Connes at the level of trees and Feynman 
diagrams [121 1 1 3] sparked a reformulation of many quantum field constructions (renormalization [141 
HS1 [TS] , Wick's theorem [T7] , quantization fTTJ [TB] , structure of Green functions [HH HU] , gauge theory 
[21]). The result of these efforts is a reasonably complete presentation of QFT in terms of Hopf algebraic 
concepts [35]. Hopf algebras, which are powerful tools to solve combinatorial problems 23J, could be 
expected to help also in the presence of a general state. Indeed, the use of Hopf algebraic methods 
resulted in the determination of the equation of motion of the Green functions |24j and the description 
of the relations between general and connected Green functions in the presence of a general state [2"2"] . 

The relation between connected and one-particle irreducible (1PI) Green functions, which is the main 
topic of the present article, is more difficult to understand and depends on the very definition of when a 
diagram is irreducible. Here, we show that a rather natural definition leads to a complete description of 
connected Green functions in terms of 1PI Green functions. 

We should point out that our main long-term interest is in the study of the electronic structure 
of highly correlated materials by means of Green functions, also in the framework of classical (non- 
relativistic) quantum chemistry. However, we restrict here our study of one-particle-irreducibility to the 
particular case of local potentials. This covers the potentials of quantum field theories -including QED, 
which is meaningful for our long-term purposes-, but not the Coulomb potential. There are several ways to 
remedy this problem and extend the constructions in the present article to time-dependent perturbation 
theory in many-body physics. For example, as the referee pointed out, our results may be extended to 
the case of general interactions by means of the notion of block truncation as introduced in [5_ . However, 
for the sake of simplicity -and since the combinatorics of one-particle-irreducibility is already intricate 
enough for local potentials, we decided to postpone the study of electronic systems (and of non-scalar 
fields) to further work. 

The paper starts with a short introduction to Green functions in the presence of initial correlations 
and to Hopf algebras, followed by the definition of quantum field forms and their convolution logarithm. 
Then, the relation between forms and connected forms is explicited, providing the classical relation 
between general Green functions and connected Green functions. To discuss 1PI functions, we need to 
generalize a recent work by Mestre and Oeckl jT^J [2D] . Then, a rather natural definition of 1PI functions 
will be proposed and the Mestre-Oeckl approach will be used to write a connected Green function in 
terms of these 1PI functions. In the process, universal properties of symmetric functions with respect to 
Hopf algebra derivations are put forward. 
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2 Green functions with initial correlations 



In the present section, we fix the notation and survey briefly the definition of Green functions, emphasizing 
the role of the initial state. 

2.1 Field operators 

The quantum fields are operator- valued distributions acting on a Fock space |25j . Here, we describe the 
construction of Fock space, creation and annihilation operators and the corresponding quantum fields. 
We start from a self-adjoint operator h acting on a Hilbcrt space Ti and, for notational convenience, we 
assume that h has a pure point spectrum, so that there is an orthonormal basis |e.;) (with i 6 /) of Tt 
consisting of eigenvectors of h. In many applications, the Hilbert space TL is a function space and the 
eigenvectors (written then preferably with the functional notation <j) n ) are functions of r (where r is a 
point in three-dimensional space). 

The tensor product of Hilbert spaces is well defined (see Ref. |26j p. 49) and the symmetric Fock space 
over 7i is the Hilbert space S(H) — ©jv=o S N CH), where S N (H) is the iV-fold symmetric tensor product 
of TL. An orthogonal basis of the vector space S N (H) is provided by the set of vectors 



where a runs over the permutations of N elements and where runs over the subset of I 

such that i\ < • • • < i^. In this formula, the symbol V denotes the symmetric product and \j\TN~\ is a 
normalization factor. 

In many-body theory, S N (H) is called the iV-particle space of the system and its elements are the N- 
particle states. In particular, S°(Tl) is a one-dimensional vector space denoted by CI in the mathematical 
literature. In many-body physics and quantum field theory, this unit 1 of the tensor product is denoted 
by |0), this is the vacuum of the theory (i.e. the state without a particle). 

The creation operator a^ n is defined as the linear map from S(H) to itself such that, for any basis 
vector \u) of S(l-t), a^|u) = c n {u) \e n ) V where c n {u) is a normalization factor (see for example [27 ). 
It is called a creation operator because it maps S N (TC) to S N+1 (H): it adds a new particle to a TV- 
particle state. Its adjoint a n is called an annihilation operator. The normalization factor ensures that 
the commutation relation a m ° o) n — a^ n o a m = S nm holds, where o denotes the composition of operators. 

In functional notation, the corresponding quantum field is the (self-adjoint) operator-valued distribu- 
tion on the three-dimensional space: 



This formalism is used to describe scalar particles or photons (up to an additional vector index in that 
case; recall that a self-adjoint field operator describes a neutral particle, charged scalar or fermion field 
operators are not self-adjoint). We remind that we focus in the present article on scalar particles (without 
a charge). 

2.2 Adiabatic limit 

The adiabatic limit is a very general way of solving the Schrodingcr equation for a system described by 
the Hamiltonian H = Hq + V where the eigenstates of Hq are known but not those of H . The general 
idea behind the technique is that, for a particle evolving in a potential V , the effect of the potential on the 
motion can be treated (or can be expected to be treated) perturbatively. As we mentioned, besides the 
perturbative expansions of QFT, the same general idea shows up in the perturbation series of statistical 
physics [5]. 

As far as adiabatic limits are concerned, the basic idea is quite simple. We define a time-dependent 
Hamiltonian H{t) = Hq + e~ e ^V. When e is small, the interaction H(t) is very slowly switched on from 
t = — oo where H(—oo) = Hq to t = Q where H(Q) = H. It is hoped that, if e is small enough, then an 
eigenstate of H is transformed into an eigenstate of H. 




(1) 



X^„(r)a„ + <(r)a; 
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To implement this picture, the time-dependent Schrodingcr equation id\$l s{t)) / dt — H(i)\^fs(t)) 
is solved. However, the solution \^>s(t)) is not convenient because it has no limit when t — > — oo. 
Therefore, we define |#(t)) = e iHot \^ s (t)) that satisfies id\V(t))/dt = H int (t)\^ {t)) , with H int (t) = 
e iH ty e -iH t e -e\t\_ j\j ow £f int (— oo) = and oo)) makes sense. Using Hi nt , we can start from the 
ground state |$o) of Hq and solve the time-dependent Schrodinger equation with the boundary condition 
|$(— oo)) = |$o)- When no eigenvalue crossing takes place, |$o) should be transformed into the ground 
state |*(0)) of H. 

Instead of calculating directly ^(t)) it is convenient to define the unitary operator U{t) as the solution 
of idU(t)/dt = H int (t)U(t), with the boundary condition {/(-oo) = 1. Thus, |#(*)) = U(t)\$ ). Note 
that U(t) depends on e, as H- lni {t). But is lim e ^o [/(0)|$o) an eigenstate of HI It would if the limit 
existed, but it does not. However, Gell-Mann and Low [25] discovered in 1951 that 

I.CL) = lim "™*> 



($„|£7(0)|$o) 



exists and is an eigenstate of iJ. A mathematical proof of this fact for reasonable Hamiltonians came 
much later [29]. Notice that the above scheme works when the ground state of Hq is non degenerate. 
When it is degenerate, the problem is more subtle and the limit e — > holds only when |$o) is well chosen 

[mini]. 



2.3 Green functions 

We now come to the heart of QFT: the calculation of Green functions (or moment functions). Green 
functions are important because they allow for the calculation of practically all relevant physical ob- 
servables: energy, charge density, transport coefficients, current density, dielectric constants, etc. In 
particular, they show up naturally in the perturbative expansions arising from adiabatic limits. If we 
could calculate Green functions exactly, we would know all interesting properties of matter. Of course, as 
far as many-body theory is concerned, we cannot calculate exact Green functions for realistic materials, 
but non-perturbative approximations are now used with great success |32j . 

When the dynamics of the particles is described by a one-body Hamiltonian Hq, the n-point Green 
function for scalar particles is defined by 

G° n ( Xl ,...,s n ) = ($ \T(tp(x 1 )...<p(x n ))\$ ) ) 

where x — (t, r), T is the time-ordering operator and ip(x) is related to v 5 s( r ) by 

<p(x) = e' ff »Vs(r)e-' ffot = ]T e" ie » Vn(r)a„ + e <e »V;(r)a+ , 

n£l 

where the ^ n (r) are eigenvectors of Hq with associated eigenvalues e„. The time-ordering operator orders 
the quantum fields ip(x\), . . . , tp(x n ) so that the field ip(xi) is on the left of f(xj) if U is greater (i.e. 
later) than tj. For example Tup(xi)(p(x2)) — <p(x\)ip{x2) if t\ > and T {ip{x\)y>(x-i)) — (p(x2)<p{x\) if 
ti < t 2 . 

When the dynamics of the particles is described by a Hamiltonian H = Hq + V, where Ho is one-body, 
the expression for the Green function becomes 12"?] 

($ \T(ip(x 1 )...<p(x n )e- i f H ™M dt )\Q ) 
n(2;i '---' SnJ ($ |T(e-/^(*) dt )|$ ) 

where H int (t) = e iH ty e -iH t e -e\t\_ Tllis g enera jj zes to 

non-scalar particles and, in the example of the 

non-relativistic electrons, 

H int (t) = e- e l*li /'drdrV t (t,r)V t (t,r')V;e(r-rOV^r / )V'(t,r). 

We are now ready to enter correlated systems. Assuming that the initial state |$o) is the free field 
vacuum |0) implies the classical expansion of Green functions in terms of Feynman propagators and, 
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ultimately, of (usual) Fcynman diagrams. In many-body theory the use of this decomposition of Green 
function into sums of (usual) Feynman diagrams is restricted to very specific states |$o) called quasi- free 
states. For the other states (or, equivalently, for the study of general functional measures in statistical 
physics) , the structure of Green functions is more complex. Let us give a simple example. We can define 
the quantity D4(xi,X2,%3,%4) by 

Gl(xi,x 2 ,x a ,Xi) = G 2 (xi 1 x 2 )G 2 {x 3l X4) + G 2 (xi, x 3 )G 2 {x 2l x 4 ) + G%(xi, x 4 )G 2 (x 2 , x 3 ) 
+D i (x u x 2 ,xz,Xi). 

When the initial state is the vacuum or a quasi-free state, the term D4 is zero. For a general initial state, 
it is not. 

For a fermionic system, a term D4 can be defined similarly. It is absent when the ground state of Hq 
can be written as a Slater determinant. It is present when the ground state of Hq is degenerate, as in open 
shell systemtQ- In that case, |$o) can be written as a linear combination of Slater determinants and D4 
describes the correlation between these determinants. The presence of several Slater determinants in the 
intial state is rather catastrophic for many-body theory. Yaris and Taylor summarized the situation [36] : 
"The inability to handle open-shell systems is a ubiquitous problem in many-body theory. It basically 
arises when one cannot find a single-determinant unperturbed ground state which connects to the exact 
ground state when the residual interaction is adiabatically switched on. When this situation holds, one 
cannot properly define occupied and unoccupied single-particle states, Wick's theorem does not hold, and 
Dyson equations, Bethe-Salpeter equations, etc. do not exist." To this list one can add that ordinary 
Feynman diagrams and the Gell-Mann and Low formula are lost. In other words, most of the tools of 
quantum field theory break down. Since the seminal work by Bloch and Horowitz in 1958 [37], many 
works were devoted to the rebuilding of these tools. Morita discovered a modified version of the Gell- 
Mann and Low theorem |38J . Fujita defined generalized Feynman diagrams [39] . Hall derived a Dyson 
equation for degenerate systems [4] . Since then, progress has been quite slow because of the combinatorial 
complexity of the problem. 

To illustrate this complexity, we first describe the generalized Feynman diagrams introduced by Fu- 
jita [SI], Hall [1] and Djah et al. [S]. For bosonic and fermionic systems, D4 can be thought of as a sort 
of 4-point Feynman propagator, as D 2 (x, y) = G 2 (x, y) is the 2-point Feynman propagator. We shall see 
that D4 plays the role of a cumulant, as in the decomposition of a distribution function. Higher order 
Green functions G 2n give rise to higher order propagators D 2n and the precise relation between them will 
be described in the following. In standard quantum field theory, the Green function of the interacting 
system can be written by adding all possible Feynman diagrams involving the two-point propagator D 2 . 
When the initial state is not quasi-free, the Green function is written as the sum of all possible Feynman 
diagrams involving 2-point, 4-point, and 2n-point propagators for arbitrary n. An example will be given 
in figure [2] of this paper. 

3 Hopf algebra 

We do not provide the general definition of a Hopf algebra [30], but we consider only the special case 
of the symmetric Hopf algebra S(V) = Q)S n (V) = (J) V®" /S n , where V is a complex vector space 

n n 

and where S n , the symmetric group of order n, acts by permutation on the components of the tensor 
power V® n , The commutative product of S(V) is denoted by concatenation. The counit is the linear 
map £ : S(V) -> C defined by e(l) = 1, e(u) = if u e S n (V) with n > 0. The coproduct is the 
linear map A : S(V) -> S(V) ® S(V) determined by Al = 1 ® 1, Aa = 1 ® a + a ® 1 for a e V 
and A(uv) = (Au)(Av), for u and v in S(V). We employ the strenghtened Sweedler notation for the 
coproduct [41] : Au = u (1) ® u {2) - Recall that there is an implicit summation in the notation, which does 
not lead to ambiguities when handled correctly: the right hand side should be understood not as the mere 
tensor product of two elements in S(V) but as a sum of such elements (so that e.g. a (1) <g>a (2) = l<g>a+a<8)l 
for a G V). More generally, an expression such as u (1) v {1) u (2) v (2) , which stands for A(u)A(v) — A(uv), 

1 For degenerate systems additional complications come from the fact that time-ordered products must be defined over 
a closed time path I . i 4 1 and the initial state must be chosen carefully. However, this does not change the combinatorial 
aspects of the problem. See ref. |35] for more details. 
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contains an implicit double summation and should be understood as: (tt (1) ®U( 3 ))(^(i) 8>«(2)), and similarly 
for expressions of higher orders. 

The iterated coproducts A k are defined by A = id, A 1 = A and A k+1 = (id® fe ® A)A fe . Their action 
on an element u of S(V) is denoted by A k u — <8 ■ ■ '®U(fc+i)- F° r an y u S ^(^O) the reduced coproduct 
is the map A : 5(F) — > 5(V) (8 S(V) such that An = Am— 1®u — The iterated reduced coproducts 

A fe are defined by A = id, A 1 = A and A fe+1 = (id® fc ® A)A fc . Their action on an clement u of S(V) 
is denoted by A k u = ® • ■ • <8 «(/c+i). The coproduct and the reduced coproduct are cocommutative, 
that is: 

A(u) = u ( i) <8 u (2) = u (2) (8 U(D, = "a) <8 Ugj = u ( 2) <8 Uq). 

The coproduct is an algebra morphism, but the reduced coproduct is not. Its relation with the product 
is described by the following simple and useful lemma. 

Lemma 3.1. If a £ V and u £ S(V), then 

A (an) = a<8u + u<8a + aw a) 8> u {1) + u a) <8 au {1) , 

and, for k > 1, 

A k (au) 

More explicitly, for k > 0, 
A k (au) 



where the terms i = 1 and z 
(l® fc <8 a)A k u in the second term. 

For an arbitrary u £ S n (V), n > and v £ S(V), we also have: 

A k (uv) = u w v w ® ••• ® u ih) v m 

= Y v (1) ®---(g>u { ^ ) v (tl) (E>---®u ( ^ ) v (tp) ®---®v (k) . 

l<p<k l<ii<— <i p <fe 

4 Green functions for quasi-free states 

Let V be the vector space generated by the symbols tp(x), where x runs over points of Mr. In physical 
terms, tp should be thought of as a free bosonic field operator, that is, as an operator-valued distribution 
(think of the quantum fields 4>s (r)). Our forthcoming developments can be adapted easily to fermionic 
systems, the adaptation amounting mathematically to replacing the symmetric algebra S(V) by the 
exterior (or Grassmann) algebra A(V"), see [TTj . 

Defining a time-ordered product of fields at the same point gives rise to major difficulties and is the 
subject of renormalization 42J. Here, we take advantage of the fact that the combinatorics of Green 
functions is in many respects a self-contained topic and leave aside these questions (renormalization, 
operator product expansion). We will therefore treat powers of fields such as tp 4 (x) as formal expressions, 
that is as monomials belonging to the symmetric Hopf algebra S(V) D S' 4 (V r ). Note that tp (x) = 1 is 
the unit of the algebra S{V)E 

2 In another paper |22] , an algebra different from S(V) was used, where tp°(x) was not the unit of the algebra, in order to 
obtain some desirable coalgebraic properties. That alternative construction considers the field products ip n (x) as the basis 
of a Hopf algebraic fibre at x. However such a point of view is not required in the present paper. 



a ® A fc ~u + aU(T) 8) A fe 1 u ( 2 ) + u a) 8) A fe 1 (au (2) ). 



fc+i 



^ u a) <8 • • • <8 W(i-i) <8 a (8 <8 • • • (8 U(jg 

i=l 
k+1 

+ u <±> ® " • • ® u (i-i) ® au(i) (8 <8 • • • 8) U(fe+i) , 

i=l 

k + 1 are a eg) A fc u and A fc w 8) a in i/ie /irsi sum and (a 8) l® fc )A fe -u and 
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4.1 Convolution 



In this section we survey some Hopf algebraic concepts that provide a startling simplification of the 
decomposition of the expectation value of time-ordered product in terms of Feynman diagrams. We 
define a form as a linear map from S(V) to C. A unital form is a form p such that p{\) = 1. In our 
context, that is when (p(x) is the quantum field of QFT or many-body theory, unital forms are defined 
from states of Ho: if |$) is a normalized state and u £ S(V), then p(u) = ($|T(u)|$) is a unital form 
because it is obviously linear and p(\) — ($|1|$) = 1. The unital form corresponding to the vacuum is 
denoted by po, so that po{u) — (0|T(m)|0). 

To express po{u) in Hopf algebraic terms, we first need a few definitions. The convolution product 
of two forms p and a is the form p * a defined by (p * f = p{u {l) )a(u i:2) ). Notice that, because of 
the commutativity and cocommutativity of S(V), a * p = p * a. The space of forms equipped with the 
convolution product is a commutative group, denoted by S, whose unit is the counit s. 

The n-th convolution power of a form p is the form p* n defined recursively by p*° — e, p* 1 — p and 
p*{n+i) _ p*n ^ p rpj^ convo i u n on exponential of a form p is the form e* p defined by 

n=0 

The convolution logarithm log* p of the form p is the form defined by 

iog> = Y,^J—( p - e y n . 

n=l 

Note that, if p is a unital form, log* p satisfies log* p(l) = 0. A form a such that <r(l) = is called an 
infinitesimal form because it is the logarithm of a unital form. In the present paper, the convolution 
exponential is always applied to infinitesimal forms. Note that, if a = log* p, then e* CT = p. In other words, 
convolution exponential and convolution logarithm are inverse functions of each other. At last, note that, 
if a and [3 are two unital forms with convolution logarithms a and &, then a * /3 = e* a * e*'' = e*( a+b \ 

4.2 Expansion in Feynman diagrams 

In standard quantum field theory, Wick's theorem states that, if u = ip kl (x\)...ip kn (x n ), < Q|T(w)|0 > is 
calculated as the sum of all pairings of k\ times the point x\, . . . , k n times the point x n . A pairing is 
the choice of a pair of different points represented graphically as a line and analytically as a Feynman 
propagator. Graphically, po(u) — (0|T(w)|0) is therefore represented by the sum of all the graphs with n 
vertices labelled by xi, . . . ,x n such that hi edges are incident to the vertex labelled by x i: for i = 1, . . . , n. 
Each graph is weighted by a proper combinatorial factor. 

To express pa(u) in Hopf algebraic terms, we define the infinitesimal form r by 

T(ip(xi)ip(x 2 )) := D F (x 2 - xi) if xi ^ x 2 , 

and 

T(cp(xi)...ip(x n )) := 0, if n =/= 2 or n = 2 and x\ = x 2 - 
The form r is called the Feynman form. The function Dp is definecj^l by 

Df(x)= [ ^- 4 - 2 V— ^ (P ' X) - 

w J (2tt) 4 p 2 - m 2 + le 

We can now restate Wick's theorem algebraically: 

Theorem 4.1. (Jb \ [j] [771 \22$ ) The unital form po is the convolution exponential of the Feynman form : 

p = e . 

3 When xi and X2 are separated by a light-like interval, the definition of r does not make sense and Dp should be 
replaced by a smooth rcgularization. We do not enter into these details here because we consider only the combinatorial 
aspects of the problem. 
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This theorem extends to the case of quasi- free states [43] , the only change is that r is now defined by 

T(ip(xi)ip(x 2 )) := D 2 (xi,x 2 ) if xx ^ x 2 , 

and 

T(ip(xi)...ip(x n )) :— 0, if n 7^ 2 or n = 2 and xi = 22, 
where D 2 (xi,x 2 ) = ($|T(<p(xi)<p(a;2)) |$). 

5 Green functions for general states 

Most often the relevant object to deal with in perturbative expansions is actually not the unital form p 
built from the ground state |$) (or, abstractly, the group S) but its convolution logarithm r (resp. the 
corresponding commutative Lie algebra £). The infinitesimal form r is called the cumulant form. As we 
shall see in section [621 this is exactly what we need to calculate the Green functions in the presence of 
initial correlations. 

The theorem (|4.ip of the previous section trivially generalizes to arbitrary states and unital forms: 
Theorem 5.1. f|3 [7| The unital form p is the convolution exponential of its cumulant form : 

p= e . 

Although our approach is not the usual one, writing p as e* r is in fact quite common in physics. 
The notion of a cumulant form is related to the cumulant expansion, and expresses the generalized Wick 
theorem used for solving the Anderson model [2]. Moreover, it is a way to isolate the singularities of the 
forms because a natural property of a quantum field is that r{(p{x\) . . . (p(x n )) is a smooth function of 
xi, . . . ,x n , except possibly for n — 2 (see [55], as well as [J5] for a related result in many-body theory). 

Finally, an observation that will prove essential in our forthcoming developments: all our previous 
reasonings suggest that a unital form p = e* r should be dealt with by means of generalized propagators 
in the same way as vacuum expectations of time-ordered products of free fields arc dealt with by means 
of 2-points Green functions and Feynman propagators in the usual picture of QFT. However, whereas 
the Feynman propagator, which is associated to the unique non trivial component of r (recall that r = 
on S l (V) for any i 7^ 2), is described graphically by a line linking two vertices, we may have now 
r(<p(xi) . . . ip(x n )) ^ with n 7^ 2. Accordingly, we shall represent graphically the "n-point propagator" 
D n (xi, ... , x„) = r(tp(xx) . . . </?(x„)j3 by a white dot with n edges linked to the n vertices x\ , . . . , x n , as 
shown in figure[T](a similar convention was used by Djah et al. [5]). 

D 3 (x,y,z) = 

x y 

Figure 1: The generalized propagator D$(x, y, z) — rUpix), (p{y), <fi{z)) 

As we already mentioned, we can also consider p = e* r as a generalization of Wick's theorem when 
the latter is stated algebraically. The same observation holds for graphical statements of the theorem: 
we saw that, in standard quantum field theory, Wick's theorem states that, if u = ip kl (xi)...ip kn (x n ), 
(0|T(zt)|0) is calculated as the sum of all pairings of k\ times the point x\, . . . , k n times the point x n . In 
the many-body context e* T (u) is replaced by e* r (u). This amounts to say that we write e* r (u) as the sum 
of all ways to partition the multiset made of k\ times point x\, . . . , k n times point x n into sub-multisets 
of any multiplicity (i.e. not only pairs and not only different points). See figure [2] for an example. The 
n-point propagators are then a convenient way to represent these sub-multisets. 

To conclude this section, we state three easy but important lemmas 

4 The definition of D n generalizes the definition of D4 in section [2 . 31 -this should be clear from our forthcoming develop- 
ments. 
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Lemma 5.1. If a G V and E = e Xa , then for any form p with logarithm r, we have p(E) — e r ( E \ 
Lemma 5.2. If a G V and u G S(V), then, for any linear map r : S(V) — ► C such that r(l) = 0, 

e* r (au) = y"V(aM (1) )e* r (M (2) ). 

More generally, 
Lemma 5.3. For any u G kere and any v in S(V), 



e* r (uv) = ^yr(u a) v w )...r(u (k) V( k) )e* r (V( k+l) ). (2) 
fe=i 

Proof. The first lemma is a simple consequence of the fact that E is group-like (that is, A(E) = E E). 
The second lemma was shown in 22 , it follows from the cocommutativity of the coproduct and the fact 
that a is a primitive element (that is, A(a) = 1 (g) a + a ® 1). The third lemma follows from the last 
identity in Lemma l3.ll from the properties of the binomial coefficients, and from the cocommutativity of 
the coproduct. □ 

The first lemma is often used with a = J j(x)ip(x)dx (up to a suitable extension of the definition of 
V when the function j(x) has not a discrete support). In that case, it relates the generating function of 
the moments of p to that of the moments of r. The second and third lemmas provide powerful tools for 
the recursive proof of the properties of e* r . Notice in particular that, using the last Lemma with v = 1: 

fe=l 21 H Hfc=n 

where, for p = 1, . . . , k, i p > 0. A formula with less terms can be given using the Faa di Bruno coefficients: 

pCa «) = V ^(« 1 )° 1 -r(o")°- (4) 



ai!(l!) Ql Q!2!( 2! ) a2 ■ ■■OL n \{nX) c 

where (ai, . . . , a n ) are nonnegative integers such that ion = n - For the partition represented by a, n 
is cut into k — J2i a i parts. For example, p(a) = r(a), p{a 2 ) = r(a 2 )+r(a) 2 , p{a 3 ) = r(a 3 ) + 3r(a)r(a 2 ) + 
r(a) 3 . 



6 Connected forms 

In quantum field theory, an important simplification comes from the fact that all physical quantities can 
be expressed in terms of connected diagrams. We first define the notion of a connected form by analogy 
with that of a connected diagram. A monomial of S(V) can always be written u — ip ni (xi) . . . (p nk (xk), 
where all points Xi are distinct. Now, for any form p with convolution logarithm r, we use eq. (HJ) to 
expand p(u) in terms of r. 

Proposition 6.1. We have 

^)= e *» = e^ e 

• „i + ... +n i = „.i=l »•••• «• 
i=l...k 

x T (<p< (asi) . . . <p n *(x k ))...r(tp n 'i ( Xl ) . . . <p<(x k )), 

where, for i — 1, . . . , k, the sum is over all the l-tuples of nonnegative integers (n\, . . . , n\) such that 
nj + h n\ = rii . 
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Although a straightforward application of the Hopf algebra formalism, the result is important since 
it allows us to compute the multiplicity of a graph -or symmetry factor- in the Feynman diagrammatic 
perturbative expansion of amplitudes. We refer for example to the expansion of the connected Green 
functions for (p 3 theory with an arbitrary ground state in the present section of the article. 

Let us consider a term t := r(</? n i (xi) . . . ip nk (xk)) ■ ■ ■ r(ip n i (x±) . . . ip nk (xk)) of e* r (u). We say that 
xi =j Xj,l < i,j < k if there exists m < I with ri^nj 1 ^ 0. The transitive closure =t of the bi- 
nary relation = t dehnes the connectedness of t: t is said to have n connected components if there 
are n equivalence classes associated to the equivalence relation = t . The connected component of Xi 
in t is defined similarly as the product of all the r((p ni (xi) ... </?"<= (xk)) with n" 1 ^ for at least 
one coefficient j with x\ =% Xj. When n = 1 (resp. n ^ 1), we also say that the term t is con- 
nected (resp. disconnected). Let us take a simple example. For u = (p(x)tp 2 (y), we have p(u) — 
r(ip(x))r(tp 2 (y)) + r{ip{x))r{ip{y)) 2 + 2r((p(x)tp(y))r(ip(y)) + r(ip(x)(p 2 (y)), where the first two terms are 
disconnected (they actually have two connected components). The connected components of y in the 
four terms are respectively r(tp 2 (y)), r(ip(y)) 2 , r((p(x)tp(y))r(ip(y)) and r((p(x)ip 2 (y)). The definition of 
connected form is actually best formulated in algebraic terms: this is the purpose of the next section. 

6.1 Another coproduct on S(V) 

As we have just seen, a pedestrian definition of connectedness makes an essential use of the fact that 
some points Xi are equal or distinct. In order to reflect this distinction, we define a new coproduct, the 
disconnecting coproduct 5 : S(V) — + S(V) (S> S(V). 

So we write a monomial of S(V) as u — p ni (xi) . . . ip nk (xk) 1 where all points Xi are distinct, and 
we define the coproduct of u as follows: 5ip n (x) = 1 ® <p n (x) + <p n (x) <8 1 if A; = 1, and Su = 
5(ip ni (xi))5(ip n2 {x2)) ■ ■ ■ S(ip 7lk (xj;)) if k > 1. Notice that this coproduct is coassociative and cocom- 
mutative but is not an algebra morphism, because 8(<p 2 (x)) ^ (8(p{x))) 2 . Since 5 is coassociative and 
cocommutative, we may still define an associative, commutative and unital product *, the disconnecting 
convolution product, on Lin(S(V), C): 

V(/, g) £ Lin(S(V), C), fig := tt o (/ g) g) o 5, 

where 7r denotes the product: 7r o (/ ® g) (u ® v) — f(u)g(v). The unit of * is the same as the unit of * 
(the projection map e from S(V) to C C S(V)). To distinguish between the two products * and S, we 
write the operations involving i with a superscript S: for example, we write log*, and so on. 

The relation between 5 and A is investigated in 22J. The reduced coproduct 6 and the iterated 
coproduct 8 k are defined as in section[31 The enhanced Sweedler notation for the disconnecting coproduct 
is Su = u {1} ® u {2) and 8u = u {L} ® U{2>- 

The new coproduct S enables us to give an algebraic definition of the connected form p c corresponding 
to the unital form p: 

p c = log*(p) 

that is, 

Vu e S(V), Pc (u) = ~^—(P ~ £ X"m)- • -(P ~ e)(« <n} ). 

n— 1 

This identity can be understood as a linked cluster theorem for forms. The same problem was addressed 
in [SI Sect. 6] with another approach: in that article, the authors introduce two notions of truncated 
moment functions, indicated respectively by an exponent T and (T). This corresponds roughly, at the 
Hopf algebraic level, to our distinction between the two coproducts A and 5. We thank the referee for 
pointing out to us this point and refer to [S] for further insights on truncated moment functions. 

Pedantically, the set of connected forms is defined as the image of the group of unital forms under 
the map log*. The two sets are in bijective correspondence and, reciprocally, we can express any unital 
form p in terms of p c by 

p — e Pc 

or: 

Vti S S(V),p(u) = e{u) + ^2 -|PcKi}) ■ ■ ■ Pc{u {n} ). 

n=l 
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For example, p c (<p n (x)) = p(<p n (x)) for n > and p c {f m (x)ip n (y)) = p(ip m (x)ip n (yj) ~ p((p m (x))p(<p n (y)) 
for to > and n > 0. The connected form p c is an infinitesimal form (that is, p c (l) = 0). For u as above, 
p c (u) is defined as the sum of the connected terms of p(u). This terminology is due to the fact that we 
can define Feynman diagrams to represent p[u), and that p c (u) is obtained by summing the connected 
Feynman diagrams present in p(u). 

Note that the relation p = e* Pc is the analogue of the relation Z = e w between the partition function 
and the free energy. As we prove now, the latter is a consequence of the former. The Hamiltonian density 
of a quantum field theory of local interactions is of the form H(x). Therefore, SH(x) — H(x)<g)l+l&)H(x) 
is primitive and, by lemma 5.1, Z = p(S) = e , where W = p c (S) (see 22 for a detailed proof). This 
is a generalization of the standard relation to the case of a general initial state. As a matter of fact, 
the coproduct S was precisely defined for the Hamiltonian to be a primitive element. This ensures the 
standard relation between the partition function and the free energy (compare to [5]). 



6.2 Example of the tp 3 theory 

As we saw in section 12.31 the two-point Green function for a system described by the interaction Hamil- 
tonian density u = (p 3 (x) is given by the expression: 

_ (0\T(y(xM y )e- iu )\0} 
[ ,V> (0|T(e-™)|0) 

We recall that the denominator cancels the divergence of the adiabatic switching of the interaction. In 
graphical terms, the denominator (0|T(e _m )|0) is the sum of all the vacuum Feynman diagrams (i.e. the 
diagrams that are linked neither to x nor to y) . Another way to obtain a convergent expression is to use 
the connected Green function G c (x,y) which is the sum of all the connected diagrams in G(x,y). 

For a general form, the factorization of the adiabatic divergence is more complex [331 E01 \M\ and it 
holds only for specific initial states |$) and by using a modified definition of the time-ordered product 
called the Keldysh approach [34], which will not be considered here. For notational convenience, we do 
not write the denominator in the definition of the Green functions for a general form and we put 

G(x,y)=p{<p{x)<p(y)e- i (5) 
The connected Green function is defined as 

G c (x,y) = Pc (<p(xMy)e~ l f-~ H ^). 
The term Hi nt (t) can usually be written J dxP(x), where P{x) is a polynomial in (p(x). Therefore, 

(—i) n f 

G(x,y) = p(tp(x)ip(y)) + ^ — / dxx . . . dx„p(ip(x)ip(y)P(x 1 ) . . . P(x n )) . 

n=i n ' J 

For example, in the ip 3 theory, we have P(x) = i/? 3 (a;)/3! and the first terms of the total Green function 

are 

G(x, y) = p(ip(x)p(y)) - ^J dx 1 p(ip(x)<f(y)ip 3 (xi)) - / dxidx 2 p((p(x)(p(y)if 3 ( Xl )if 3 (x 2 )) + ■■■ 

For notational convenience, we assume that D n (xi, . . . ,x n ) — r{tp{x\) . . . ip{x n )) — if n is odd. The 
connected Green function G c (x, y) is obtained by keeping the connected terms of the total Green function. 
In standard quantum field theory, the expansion to the second order gives us 

G c (x,y) = D 2 (x,y) - - / dz dw D 2 (x, z) D 2 (y,w) D 2 (z,w) 2 + . . . (6) 
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For a general form, the expansion to the second order gives a more complex result: 
G c {x,y) = D 2 (x,y)-J dz dw (^D s (x,y,z,z,z,w,w,w) + ^D 2 (x,z) D 6 (y,z,z,w,w,w) 
+ -^D 6 (x,z,z,w,w,w) D 2 (y,z) + ^Da{x,y,z,w,w,w) D 2 (z,z) 
+ ^D e (x,y,z,z,w,w) D 2 (z,w) + — D 4 (x, y, z, z) D 4 (z,w,w,w) 

+—D 4 (x, y, z, w) D 4 (z, z, w, w) + -D 4 (x, z, z, w) D 4 (y, z, w, w) 
8 4 

+ -^D 2 (x,z) D 2 {y,z) D 4 (z,w,w,w) + -D 2 (x,z) D 2 (y,w) D 4 (z, z,w,w) 

+ -D 2 (x, z) D 4 (y, z, w, w) D 2 (z, w) + -D 2 (x, z) D 4 (y, z, z, w) D 2 (w, w) 

+ ^D 4 (x,z,w,w) D 2 (y,z) D 2 (z,w) + ^D 4 (x,z,z,w) D 2 (y,z) D 2 (w,w) 

+ ^-D 4 (x,y,w,w) D 2 (z,z) D 2 (z,w) + -D 4 (x,y, z,w) D 2 (z,z) D 2 (w,w) 
4 o 

+-D 4 (x,y,z,w) D 2 (z,w) 2 + ^D 2 (x.z) D 2 (y,z) D 2 (z,w) D 2 (w,w) 
+^D 2 {x, z) D 2 (y, w) D 2 (z, w) 2 ^ +... (7) 
These terms can be given the diagrammatic representation of figure [21 

Gc( „) = _ + x. + a< + >'\ + 2< + .ss. + X 



.>uX + >€ + .s. + X 



Figure 2: The first few terms of G c (x, y) 

It is clear from a comparison of equations ([6|) and ([7|) that the use of a form p which does not come 
from a quasi-free state increases significantly the combinatorial complexity. As a consequence, little was 
known about the structure of the connected Green functions in the general case, especially as far as 
one-particle-irreducibility is concerned. 

7 Symmetric functions and derivations 

In [HI [20], Mestre and Oeckl have proposed a powerful Hopf algebraic tool to generate 1PI diagrams. 
The forthcoming developments of the present article aim at extending their work to correlated systems, 
and at describing in this general setting the decomposition of connected Green functions into one-particle 
irreducible Green functions. 

In order to do so, the langage of functional derivatives proves to be a very convenient framework. This 
leads naturally to the study of Hopf algebra derivations, which is the main topic of the present section. 
Namely, we describe some aspects of the algebra of symmetric functions in the Hopf algebraic setting and 
prove that this algebra has universal properties with respect to derivations acting on commutative Hopf 
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algebras. We will show later that this formalism nicely encodes the properties of the generalized n-point 
propagators associated to arbitrary states, as described above. We do not seek the outmost generality in 
our constructions, but mention that they can be easily extended to more general (e.g. noncommutative) 
situations. 



7.1 Definition of the derivations A rn . 

Let us consider, once again, the Hopf algebra U = S(V). Two products arise therefore in S(U) = S(S(V)): 
the symmetric product in S(V), denoted by juxtaposition, and the symmetric product in S(U), denoted 
from now on by V to avoid any confusion. 

To a unital form p on U = S(V), we associate its convolution logarithm r and, for f(xi) . . . ip{x m ) e 
U, we recall that D m (xi, . . . , x m ) — r(ip(x\) . . . tp(x m )), so that the generalized propagators D m are 
symmetric functions of their arguments. 

Recall that V is spanned by the symbols (p(x), where x runs over the points in M. d . We choose 
an arbitrary total order on the points in M. d , for example the lexicographical order on the c?-tuplcs of 
coordinates. The operators A m -i : U — ► S m (U) are then defined by 

A m -i{u) = D m(xi,...,x m ) f U{1} V ■■ 

1 \ - , . du, 1} 

Xl,...,X m 

These operators on U extend uniquely to derivations on S(U), still written A m _i, so that: 

A m _! : S k (U) — ► S k+m -\U) 

(recall, for completeness sake, that a derivation D satisfies the Leibniz rule, D(xy) = D{x)y + xD{y) 1 
so that the action of a derivation on S(U) follows by induction on the degrees from the knowledge of its 
action on U). Notice also that, by definition of S, the terms with the last summation vanish. 

An important point for our forthcoming developments is that the A^ commute. This follows from 
the Schwarz commutation rules for derivatives and from the definition of the disconnecting coproduct 
5. Graphically, these operations will allow us to construct inductively 1PI graphs (with a given set of 
vertices and with given multiplicities), and their algebraic properties will allow us to enumerate these 
graphs. This will be the purpose of the next section, whereas the following subsections introduce the 
abstract Hopf algebraic framework suited for this enumeration. 



d(p(x m ) 

v du {m} 
d(p(x m )' 



7.2 The Hopf algebra of symmetric functions 

Let X = {x\ . . . x n . . . } be a countable alphabet, and C[LY]] the algebra of formal power series over X. 
The group Soq — lim<S n acts on C[LY]] by permutation of the letters of X] the algebra of symmetric 

functions Sym is the subalgebra of Soo -invariant series in C[LY]]. 

The algebra Sym is (up to completion with respect to the filtration induced by the grading of sym- 
metric polynomials by their degree) a free commutative algebra over various families of generators. For 
our purposes, the most interesting ones are the families of power sums symmetric functions and complete 
symmetric functions associated respectively to the series 



P. :=£ft:=l + £ r 

kefi ieN 



and 



C :=]>>:=n r 



In view of our forthcoming computations, it is actually convenient to work with an extension of Sym, 

qSym: we write Qk for k > 1, Qo := q and Q. := Qk, where q stands for an additional free 

fceN 
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variable. The series are related by the familiar Newton- type identity: C. = e^*^ . We write S. for the 
^-series S. := e Qa C. = e Q * 

The algebra qSym carries a natural notion of grading (by the degrees of symmetric polynomials, with 
deg(Q n ) — n), but it is convenient, for our purposes, and for reasons that will become clear later, to 
introduce an extra "auxiliary" grading by considering the family of the Qfc as a family of generators of 
Sym of auxiliary degree 1. This is best explained through an example: Q0Q3Q5Q9 is of degree 17 and of 
auxiliary degree 5. The auxiliary degree is indicated with a superscript (whereas the degree is indicated 
by a subscript), so that, for example, the component of degree n and auxiliary degree k in S. is given by: 



where the sum runs over all (n + l)-tuples of integers a — (ao, ...,a n ) with cxq + • • • + a n = k and 
a\ + 2«2 + • • • + na n = n. Notice that we distinguish carefully between and the latter standing 
for the k-th power of S n . The following examples will be useful in the sequel: S° = S n ol, S§ = Qa/kl, 
S l n = Q n , Si = Q k Q - x Q x /{k-l)\, S| = Q Q 2 + Q?/2. 

Generating series are a useful tool to handle computations with the Qjj and the S„. Consider for 
example the series S.(a + b) := e( a + & )Q« ; its expansion as a series in the variables a and b yields: 

Proposition 7.1. We have, for all k,l > 0: 

'h + l 



n / j ^ra^n—m' 

m=0 



In particular 



1 

Cfc _ i n ofc-l 

n ~ h ' I n-m- 



m=0 



The Hopf algebraic properties of symmetric functions were recently exploited with great profit by 
Fauser and coll. [48, 49J. Similarly, we put a simple Hopf algebraic structure on qSym uniquely defined 
by requiring the power sums symmetric functions (i.e. Q n for n > 0) to form, together with Qo, a series 
of primitive elements or, equivalently, by requiring the series S, to be a group-like element. In other 
terms, the coproduct A on qSym is fully specified by requiring that A(Q n ) := Q n ® 1 + 1 ® Q n , n > 0. In 
particular, the coproduct is compatible with the two graduations. When expliciting this property with 
the S^, we get: 

Proposition 7.2. The coproduct o/S* is 

n k 

and its iterated coproduct is 



m=0 i=0 



k—i 
m ^ ^ n — m ' 



TiiH \-n p =n 

Note that propositions 17.11 and 17.21 still hold if the variables Q n do not commute. 



7.3 On Hopf algebra derivations 

Let H = ® H n be an arbitrary connected graded commutative Hopf algebra and A^^Ai, ...,A n , ... an 

ra£N 

arbitrary series of degree 0, 1, n, ... derivations on H. That is, for any p, the restriction of A n to 
Hp is a linear map from H p to H p+n , and A n satisfies the Leibniz rule: for any h, I in H, A n (h ■ I) = 
A n (h) ■ I + h ■ A n (l). We also assume that the A n commute, so that the A n generate a commutative 
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subalgebra T> of End(H) (for the composition of maps). Of course, we have in mind the particular 
derivations A m acting on S(U), but the following results hold in full generality. 

There is therefore, since qSym is free over the Q n , a universal algebra map from qSym to End(H) 
obtained by mapping Q n to A n . We write L, for the image of S, under this map, and L„ for the image 
of S„. Note that, for any p, L n maps H p to H p+n . Of course, the identities that hold in qSym for the 
Q n and the S n also hold in End(H) for the A n and the L^. More surprisingly however, the coalgebra 
structure of qSym reflects the action of D on H. We refer to [501 116] for similar phenomena occuring in 
the study of Lie idempotents and renormalization in perturbative QFT. 

Proposition 7.3. We have, for any X € qSym and any h, h' £ H : 

(3(X)(hti) = (3(X {1) )(h)(3(X (2) )(ti) 

where Xm <8> X( 2 ) stands, as usual, for the coproduct of X in qE>ym. 

The identity can be generalized by a straightforward recursion to compute @{X)(hi...h n ). Notice first 
that the identity in the Proposition is obvious when X is a Q n , since f3(Q n ) = A n is, by hypothesis, a 
derivation. Now, assume that for X and Y in qSym and arbitrary h, h' , 1,1' E H the above formula holds, 
that is: 

P(X)(hti) = {3{X (1) ){h)p{X (2) )(ti), p(Y)(W) = /3(Y- (1) )(0/?(Y(2))(O- 

It follows that: 

[3{Y)op{X){hh') = (3(Y)((3(X (1) )(h)f3(X (2) )(h')) 

= p(Y {1) ) o f3(X (1) )(h) 0(Y (2) ) o p(X {2) )(h') 
= f3((YX) {1) )(h)(3((YX) {2) )(h'). 

In other terms, if two elements in qSym satisfy the identity in the Proposition, their product also satisfies 
the identity. Since the Q n satisfy the identity, and since their products span qSym, the Proposition 
follows. 

Let us consider the particular example of S(U) and of the derivations A m -\. Let us write A m -\(u) = 
Mm— 1,1 • ■ ■ !im-i,mi with a enhanced Sweedler-type notation for the action of A m -\ from U to S m (U). 
We get: 

n+1 ^ n+1 

m— 1 m— lniH \-Tl m = n— m 

felH hfc m =fe-l 

where V is the symmetric product in S(U). This is a generalization of lemma 13 in [19] and of proposition 
15 in 20J where the Authors, Mestre and Oeckl, studied the case where A n — for n ^ 0, 1. Note that 
our notation is different from theirs. 



8 One-particle irreducible decompositions 

We consider now the derivation of an explicit decomposition of a connected Green function into 1PI 
Green functions for a general state. 

As we have just noticed, the map Q m >— > A m enables us to define L n acting on S(U) = S(S(V)). For 
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example Lq(u) = u, L^u) = A m (u) and 



Ll{u) 



2 £ Dl(aa) M*i)M* 9 )' 



dip(xi) dip(x 2 )dip(x 3 ) 



Ll{u) = \ D 2 (x u x 2 )D 2 (x 3 ,x i )-^- 



. v g 2 M{2 } v <9u {3} 



2 x 1 ,xt^3,x 4 ' ' d<p{x 2 )d<p{x 3 ) 9^(2:4) 

{3} 



L\{u) = \ £ J 1 (x 1 )Z? 1 (x 2 )g 2 ( g3 ,x 4 ) /) , V _^ 2 ' 



8.1 A tree interpretation 



The operator can be written as a sum over all the bipartite trees with k white vertices and n + 1 black 
vertices. This description in terms of trees is important because, in standard QFT, a connected Green 
function can also be described as a tree of 1PI Green functions -a description we want to extend to the 
case of initial correlations. To give a more precise relation between L\ and bipartite trees, we consider 
the expression for in terms of partitions a: 



Act„ 
sin 



a,, 



where the sum runs over the sequences a of nonnegative integers with ao+- ■ -+a n = k and a\+- ■ -+na n — 
n. The monomial corresponding to a given a is represented by the sum of all bipartite trees with k white 
vertices and n+1 black vertices, such that on white vertices have valency i + 1, for i = 1, . . . , n. 
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The terms of lowest degrees are 



T ± 



L\ = 




o — 

Li = 

L\ = 

Lq = ^J^o =^0) 

L\ = iMa + -A? =£s°-»-°- 

L\ = \aIA x 

L\ = A 3 =X, 

Lq = =( ^bi 

L 1 — 




= A 5 ="&. 

To calculate the value of a tree of L\: (i) Associate to each of the k + n edges a variable Xi, with 
i = 1, . . . , k+n. (ii) To each white vertex v, associate the factor D m {x,i x , . . . , Xi m ), where m is the valency 
of v and x^ , . . . , Xi m are the variables associated to the edges incident to v. (iii) There are n + 1 black 
vertices. Split u into n + 1 parts by 8 n u = u {1} ® • • • ® u { „ +1} . Number the black vertices from 1 to n + 1 
and to vertex I associate the factor 

d m u {li 



dip{x il )...dif>(x il )' 



where x^ , . . . , x% t are the variables associated to the edges incident to the black vertex number I. (iv) 
Multiply the factors corresponding to the black vertices with the product V in S(U). (v) Divide the 
resulting value by the order of the symmetry group of the tree. 



8.2 The 1PI components of forms 

The last step before we can write a connected form in terms of 1PI forms is to give a reasonable definition 
of what is the 1PI component of a form, similarly to the definition of the connected components of forms. 
Several definitions are possible. The simplest one was proposed by Hall [4] and has recently provided 
detailed structural results [3S]. Here we consider a definition which is strictly more general than Hall's 
and that leads to an interesting structure. In a graph, it is easy to describe what we mean by cutting a 
line or a set of lines; this approach leads, in classical QFT (with 2-point Feynman propagators) to the 
definition of 1PI Feynman diagrams as connected diagrams that are still connected when an arbitrary 
propagator line is cut. We propose to generalize the notion by replacing the Feynman form (that is, the 
classical case where only 2-point Feynman propagators are considered) by an arbitrary unital form. 

Our approach is rooted in the Hopf algebraic picture of QFT. Notice however that our constructions 
could be translated mutatis mutandis in the langagc of functional derivatives. For example, the derivatives 
0^?j that we have used in the definition of the operators Aj were defined as usual derivatives (in the 
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polynomial algebra over the symbols <p(x)) but could be understood alternatively as functional derivatives. 
The same observation holds for our forthcoming constructions. 

In proposition 16. 1[ for any u = ip ni (xi)...ip nk (x k ), we have expanded p(u) as a linear combination 
of terms such as r(tp n i{xi) . . . tp 11 * (x k ))---r((p ni (x\) . . . ip 71 * {x k )). Let us consider a connected term t := 
r((p n * (xi)... tp nl « (x k ))...r{(p n 'i (xi) ...ip nl « (x k )) in e* r (u). 

Definition 8.1. The connected term t is said to be one-particle-reducible if and only if, there exists 
i € {1, 1} and i p } C {1, k} such that 

1. (f n *(xx) ■ ■ ■ (p n ' k {x k ) = <p(x h ). ..<p{x ip ) 

2. furthermore, in the remaining part of t, 

r(p n l . . . v *l ^...r^ 1 ( Xl ) . . . p< X (x k ))r( V < +1 (xi) . . . p>< + " {x k ))...r{^< ( Xl ) . . . p< (x k )) 

the connected components of Xi ly Xi p are either empty or pairwise disjoint. 

A connected term that is not one-particle-reducible is said to be one-particle-irreducible (1PI)- 

For example, r(ip(xi)ip(x2)) 2 (a loop constructed out of two two-point propagators) is 1PI (in our sit- 
uation, and also in the usual picture), and so is r(ip(xi)(p(x2))r(tp(xi)ip(x3))r(ip(xi)ip(x2)'p(x3)), whereas 
r((p(x 1 )(p(x 2 ))r(ip(xi)p(x 3 ))r((p(x 1 )ip(x2)ip(x4,)) or r((p(xi)(p(x 2 ))r(ip(xi)(p' 2 (x 3 ))r((p(x2)p 2 (x4 : )) are not. 

Definition 8.2. The 1PI component pi of a form is the sum of all the 1PI terms in the connected 
component of p. 

Recall that forms, their connected components, and also their 1PI components are linear maps from 
U = S(V) to C. However, any such linear map I can be uniquely extended to a multiplicative map, 
still written I from S(U) to C (that is, to a character of the algebra S(U), in the algebraic terminology). 
Concretely, for u\,...,u n € U,l{u\ V • • • Vu„) :— l(ui)...l(u n ). In particular, the connected and 1PI 
components of forms can be viewed as characters of the algebra S(U), so that, for example, an expression 
such as p c o A m makes sense as the composition of a derivation of S(U) and a map from S(U) to C. 

Let us consider a simple example to illustrate these ideas, namely the structure equation linking 
connected and 1PI components of forms in the most common picture of pQFT: an interacting theory 
-say (p 3 - with Feynman diagrams built of 3-valent interaction vertices and 2-point propagators. A general 
Feynman diagram can be described as 1PI diagrams connected byn£N Feynman propagators satisfying 
the property that cutting any of these propagators makes the original diagram disconnected. Taking 
into account the symmetry factor n\ arising from the fact that these Feynman propagators can be cut 
successively in an arbitrary order results into a structure equation relating the connected and 1PI Green 
functions. In Hopf algebraic terms: 

p c = pi o e F 

where F is the derivation of S(U) associated to the Feynman propagator: for all U\, u n € U = S(V), 

with Z?2 the (quasi-)free 2-point Green function. 

In the general case, replacing Feynman propagators D2 by arbitrary propagators D n , doesn't change 
the general principles of the proof. An arbitrary Feynman diagram for an interacting theory as the ones 
considered previously in the present section can still be cut into 1PI pieces connected by a family of 
n-point propagators, n S N, in such a way that removing any of these n-point propagators splits the 
original diagram into n connected pieces. For a given n, the associated symmetry factor is p n \, where 
p n is the number of n-point propagators in the family. These observations result in a family of structure 
identities for 1PI diagrams at all orders, and an identity that should be understood as structure theorem 
for the perturbative approach to interacting theories. 
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Theorem 8.1. (Structure of connected forms) For an arbitrary connected form p c , we have 

p c = pi o L. = p! o e ^ 



£ A, 

i EN 



and 

Pi = Pc o e V »e» 

We remark that the effect of ^4o is just a shift of the fields: for instance e A ° (ip n (x)) = (<p(x) + Di(x)) n . 
For u = ip ni (xi) . . . (p np (x p ), we have A m (u) — if m > p because A m splits u into m + 1 pieces and the 
coordinates Xi of these pieces must all be different. More generally, vanishes if n > p. Because of 

the trivial effect of Aq we put Aq = and we get 

p—l niH h«p 

Pc {u) = PI { U )+J2 J2 Pl( L n( U ))- 
n=l k=l 

In particular, p c {ip n {x)) = pi(tp n (x)) and, for x ^ y, 

p c {<p n {x) V ™{y)) = Pl (<p n (x)<p m (3,)) + mnD 2 (x, y)p 1 { V n ~ 1 {x))p I ^ m - 1 {y)). 

The first equation of theorem lS.ll describes the connected Green functions in terms of 1PI Green functions. 
It is a generalization to general states of the standard QFT result and of a theorem by Mestre and Oeckl 
|20j . The second equation is new even in the QFT context: it describes the 1PI Green functions as a 
linear combination of products of connected Green functions. In QFT, 1PI Green functions are expressed 
in terms of amputated connected Green functions. Here, we do not amputate the Green functions (this is 
not allowed for a general state because parts of the Green functions belong to the kernel of the differential 
operator used in the equation of motion of the free field). 

The consequences of these identities for the QFT of interacting systems, and the fine study of con- 
nected and 1PI amplitudes are postponed to further work. 



9 Conclusion 

In this paper, we developed mathematical tools to extend the relation between connected Green functions 
and 1PI Green functions from the case of a quasi-free ground state to the case of a general state. Our 
main result is the structure of connected forms described by theorem 18.11 This work can be extended 
in two directions. On the physical side, the main structure identity can be used to derive resummation 
theorems that generalize Friedberg's formulas [51j . On the mathematical side, many of our results can 
be extended to the case of non commuting variables. 
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